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We have measured far infrared reflectance and transmittance spectra as well as Raman scattering 
spectra of a'-NaV^Os single crystals for all the principal polarizations. The temperature range 
above the phase transition temperature T c — 35 K was investigated, mainly. On the basis of this 
experimental study and of the lattice dynamics calculations we conclude that the symmetry of 
NaV2C>5 in the high temperature phase is described by the centrosymmetric D^u space group. The 
assignment of the observed phonons is given. Values of dielectric constants are obtained from the 
infrared data. Asymmetric shapes of several infrared lines as well as higher order infrared vibrational 
spectra are discussed. The crystal field energy levels of the 3d electron localized at the V 4+ site have 
been calculated in the framework of the exchange charge model using the values of effective charges 
obtained from the lattice dynamics calculations. According to the results of these calculations, the 
earlier observed broad optical bands in the region of 1 eV can be interpreted as phonon assisted 
d — d transitions. 
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I. INTRODUCTION 



The vanadate a'-NaV^Os has attracted a considerable interest recently as the second inorganic compound undergo- 
ing the spin-Peierls transition (at the highest known temperature for the spin-Peierls compounds, T c — 35 K Q). The 
spin-Peierls transition is expected to occur within a system of linear spin-1/2 Heisenberg antiferromagnetic chains 
coupled to a three-dimensional phonon field. As a result of such a coupling, magnetic atoms of the chain dimcr- 
ize and a spin gap opens One-dimensional magnetic properties of NaA^Os above 35 K follow from magnetic 
susceptibility 0, ESR || and angle-resolved photoemission 0] measurements. 

Below 35 K the lattice dimerizes, as observed by X-ray 0] and Raman J(i|,f?J scattering, infrared transmission Q 
and reflection ]9| measurements, while the magnetic susceptibility decreases isotropically, thus showing a spin gap 
formation 01. The size of the gap A = 10 meV follows from inelastic neutron scattering study of NaV^Os single 
crystals pjuj. 

The structure of NaV^Os contains double chains of edge sharing distorted VO5 pyramids running along the or- 
thorhombic 6-axis (see Fig. 1). These double chains are linked together via common corners of the pyramids to form 
the a6-layers. Na atoms lie between the layers flljjjj. The structure of NaV^Os looks like the structure of V2O5 
intercalated with sodium. In an early X-ray room temperature investigation on polycrystalline samples of NaV^Os 
Carpy and Galy [0 suggested the noncentrosymmetric space group Cl v -P2imn with two nonequivalent vanadium 
positions in the unit cell. The picture of magnetic chains of V 4+ 05 (S = 1/2) pyramids isolated by nonmagnetic 
chains of V 5+ 05 (S = 0) pyramids proposed to account for one-dimensional magnetic properties of this mixed valence 
(V 4 5+ ) compound is compatible with this space group [jlj. 

However, the recent redetermination of the structure by single crystal X-ray diffraction at room temperature was 
in favor of the centrosymmetric D^-Pmmn group with only one vanadium position in the structure | p3| , ^4[ . Though 
the topology of the structure remains essentially the same as in the previously proposed noncentrosymmetric space 
group pr[| , the possibility for charge ordering is, however, lost in the new higher symmetry group. Smolinski at 
al. |ll| and Horsch and Mack |ll| suggested a quarter-filled ladder model for NaV205, with the spins carried by V-O- 
V molecular orbitals on the rungs of the ladder. They argued that the exchange interaction along the ladder is much 
greater than that between the neighboring ladders which would explain the one-dimensional magnetic properties of 
the high temperature phase of NaV2 05. The transition at 35 K was supposed to be an ordinary spin-Peierls transition. 
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Quite recent 51 V-NMR experiment on a single-crystalline sample of NaA^Os also revealed only one vanadium position 
in the high temperature phase but pointed unambiguously to the existence of two different vanadium sites occupied 
by V 4+ and V 5+ at liquid helium temperatures |l6|. Thus, the transition at 35 K is connected with a structure and a 
charge ordering processes. Very recently, Seo and Fukuyama and Mostovoy and Khomskii jl8) proposed a zig-zag 
scheme of V 4+ -V 5+ ordering. Seo and Fukuyama argued that, as a result, two-dimensional lattice of antiferromagnctic 
dimers is formed j[7), while Mostovoy and Khomskii gave reasons in support of a system of alternating chains [fl8|| . 
Thalmeier and Fuldc |l9| have presented some theoretical reasons for the primary charge ordering which provides 
neighboring linear V 4+ and V 5+ chains with the subsequent spin-Peierls transition. Two close transitions near 35 K 
in NaV2C>5 were detected by Koppen et al. by the thermal expansion measurements p0| . 

In view of these recent works, the symmetry problem of the high temperature phase seems to be a matter of great 
urgency. Raman and infrared measurements could give an additional information to clarify the question whether the 
space group is centrosymmetric or not, because of totally different selection rules in these two cases. We reanalyzed 
our earlier infrared and Raman spectra of NaV2C>5 pl| ] and found that they are in a better agreement with the 
centrosymmetric D\\ space group than with the noncentrosymmetric C 2v one. However, in our work pl[ | we did not 
measure infrared spectra in the E||c polarization. Also, the signal-to-noise ratio of Raman spectra was rather low. In 
the present work, we reinvestigate vibrational spectra of the high temperature phase of NaV20s using different single 
crystals, including extraordinary big ones. We present far-infrared reflectivity as well as Raman-scattering spectra for 
all principal polarizations. In addition, transmittance spectra were studied. We show that our results are in a much 
better agreement with the centrosymmetric D^t group than with the noncentrosymmetric C 2v one. The assignment 
of vibrational modes is given based on a comparison with results on the previously studied V2O5 p3] and on lattice 
dynamics calculations of this work performed in the framework of the rigid ion model. 



II. EXPERIMENTAL 



Single crystals of stoichiometric a'-NaV20s used in this study were grown by a melt growth method using NaVC>3 as 
a flux |Q . Samples from different batches were used. One sample was 1.3x8x1 mm, another one was 3 x 17.3 x 1.6 mm 
along a-, b- and c-axes, respectively. For transmission measurements we have prepared four thin samples cleaved 
perpendicular to the c-axis. Their thicknesses were 110 ± 1, 45 ±5, 14 ± 1.5 and 6± 1 /im. The samples were checked 
with X-ray diffraction, magnetization and ESR measurements. They exhibited a sharp transition at about 35 K. 

Reflection and transmission measurements were performed with a BOMEM DA3.002 Fourier-transform spectrom- 
eter at nearly normal incidence of polarized infrared radiation. The following geometries of the experiment were 
used: 1) k||c, E||a and E||b; 2) k||a, E||c and E||b. Room temperature reflectance and transmittance spectra were 
measured in a spectral range 30-5000 cm -1 with a resolution 0.5-2.0 cm -1 . Using both reflectance and transmittance 
spectra, the absorption coefficient a was calculated. Low temperature (down to 6 K) transmittance spectra were 
measured with a He vapor cryostat in the spectral range 30-1000 cm -1 with a resolution 0.05-1.0 cm -1 . 

Raman spectra were excited at room temperature by the 514 nm and 488 nm lines of an Ar-ion laser in backscattering 
geometries, dispersed by a home-made triple spectrograph and recorded using a multichannel system consisting of an 
image intensifier tube with a multichannel plate and a vidicon. 



III. RESULTS 



A. Factor-group analysis 



There are two formula units and, hence, 16 atoms in the NaV2C>5 orthorhombic unit cell with lattice constants 
a = 1.1316 nm, b = 0.3611 nm, c = 0.4797 nm 1 1 1 l^|l4|]. Below, we present the results of factor-group analysis for 



both centrosymmetric D\\ (T^jlJ] and noncentrosymmetric C 2v |llj space groups. 

a. Space group D^-Pmmn The notation Pmmn refers to the standard setting, such that x\ \a, y\ \b, z\\c. It follows 
from X-ray diffraction data of Refs. jl3],[l4j that Na atoms occupy 2b positions (the corresponding fractional atomic 
coordinates are defined by the basis vectors ri( Na ) = — r2(Na) = (1/4, —1/4, zi), z\ — 0.8592) and oxygen 01 
atoms occupy 2a positions (ri(Ol) = — r2(01) = (1/4, 1/4, z 2 ), z 2 — 0.5195), both these positions having C 2v local 
symmetry. V, 02 and 03 atoms reside in different 4/ positions (i"i(A) = — r 3(A) = (xa, 1/4, za), 1*2 (A) = —r^A) = 
(1/2 - xa , 1 /4, z A ), x v = 0.40212, z v = 0.39219, x 02 = 0.57302, z 02 = 0.48769, x 03 = 0.38548, z 03 = 0.05803) with 
the local symmetry Cjf 2 . These positions yield the following irreducible representations Jj4p5| ]: 

@2v '■ T = ^9 + B 2 g + B 3g + Bi u + B 2u + B 3u 
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Cf : T = 2A g + B lg + 2B 2g + B 3g + A u + 2B lu + B 2u + 2B 3 



Multiplying the representations given above by the number of different positions of the appropriate symmetry, sum- 
marizing them and subtracting acoustic modes (B\ u + B 2u + B 3u ) , we obtain the following NaV2C>5 optical vibrational 
modes: 

T N b aV 2 oA Pmmn ) = 8A s (aa, bb, cc) + 3B lg (ab) + 
+8B 2g {ac) + 5B 3g (bc) + 3A U + 
+7Bi„(E||c) +4B 2u {E\\b) + 7B 3u (E||a) (1) 

There are 45 vibrational modes in total. A u modes being silent, 24 Raman (A g , Bi g , B 2g , B 3g ) and 18 infrared (B\ u , 
B 2ui B 3u ) active modes are expected to be found in the spectra of NaV20s, provided the crystal space group is D\\. 

b. Space group C\ v -P2\mn In their original work jnj Carpy and Galy adopted the axes setting for the P2\mn 
space group. Below, we use the standard setting for the Pmn2\ space group: x\\b, y\\c, z\\a. There are two nonequiv- 
alent V positions, five nonequivalent O positions and one Na position in this group, all of them being 2a positions 
with C| z local symmetry. In the same way as in the previous case, using tables fl25f and subtracting acoustic modes 
(Ai + B 2 + B\ ) , we find the following vibrational modes: 

r^Uo,(-P"»»i2i) = 15Ai (aa, bb, cc; E||a) + 
+8A 2 (bc) + 7B 1 (ab; E||b) + lhB 2 {ac- E||c) (2) 

There are 45 optical modes again. But in the case of this noncentrosymmctric space group all of them are Raman 
active, 37 of them are also infrared active. 



B. Infrared spectra 



Fig. 2 shows the room temperature far-infrared reflectivity spectra of NaV20s for different polarizations of the 
incident light. Experimental data are presented by open circles. Measured spectra were least-squares fitted by the 
spectra computed according to the expression 
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The classical dispersion formula for N independent damped oscillators was used: 

N 



(3) 
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For E||b and E||a polarizations the number of oscillators and initial values of parameters were taken from the 
transmittance spectra ( |2l| and the present work). The peculiarity at 1014 cm -1 in E||a polarization crossed out 
in Figs. 2 and 3 and observed also in E||b polarization for some samples depends on a particular sample and is not, 
evidently, an intrinsic property of NaX^Ov It was not taken into account in the fitting procedure. In addition to weakly 
damped phonon oscillators, an overdamped oscillator centered at about 300 cm -1 (u>j = 291 cm -1 , 7,; = 260 cm -1 , 
/, = 0.38) was introduced in E||a polarization to account for a low frequency part of a broad absorption band of a 
complex two- humped shape found in our previous study Q (see also Fig 3). We failed to model the high frequency 
hump of this band centered at about 1000 cm^ 1 by a similar oscillator and did not try to use a more complicated 
model. This results in a not so good fitting of the high frequency part of the reflectance spectrum. In E||a polarization 
the phonon at about 150 cm -1 could not be fitted well. This line is strongly asymmetric in transmittance spectra, 
obviously, due to the interaction with the underlying broad band. 

A small bump in reflection at 939 cm' 1 shown by the arrow in Fig 2 corresponds to the Fano-type resonance 
p6j well seen in the absorbance spectrum (Fig. 3). One more such a resonance becomes visible below 200 K at the 
frequency of about 91 cm" 1 (see Fig. 4 and also Ref. ||). We fitted the absorption coefficient in the vicinity of these 
two strongly asymmetric lines by the expression |2^| : 

a(u) = a B {u) + a (5) 
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where £ = (u> — w r )/7, ag(w) is a slowly varying broad band absorption (it is shown by a dashed line in the vicinity 
of the 939 cm' 1 sharp resonance in Fig. 3), oto, u> r , 7 and q being variable parameters. Such an expression describes 
various physical situations of a sharp transitions being overlapped by a broad continuum. The lineshape of a sharp 
transition is changed by an interference with a continuum and depends essentially on the strength of an interaction 
between discrete and continuum states. The parameter q being inversely proportional to the matrix element of an 
interaction, the case \q\ —00 corresponds to zero interaction and results in a normal Lorenzian resonance, \q\ = 1 
yields a dispersion-like curve, while |g| = gives an inverted Lorenzian (antiresonance) . The ratio ao/as shows what 
fraction of the continuum states interacts with a sharp excited state. The results of the fitting are displayed in the 
inset of Fig. 3 and in Fig. 4. A similar fitting should be performed for the resonance at about 150 cm" 1 but we failed 
to construct as(w) in this case. The parameters obtained from the fittings are listed in Table 1. loto an d jto stand 
for u>i and, correspondingly, 7^ of Eq. (Q) or for to r and 7 of Eq. (|5|). LO frequencies and damping constants were 
calculated as complex roots of the equation e(u) = 0. 

The left inset of Fig. 4 presents the temperature dependence of the Fano parameter q for the spectral line near 
91 cm" 1 for the temperatures higher than T c = 35 K. It should be mentioned that below 35 K, the shape of this line 
changes to an ordinary Lorenzian lineshape (see the right inset of Fig. 4 and also |^| ) . Simultaneously, a continuum 
absorption diminishes markedly in this spectral region while it practically does not change at the maximum of the 
low- frequency hump at 320 cm" 1 . 

With decreasing the temperature, besides the asymmetric resonance at 91 cm" 1 in E||a transmission, two lines at 
215 and 225 cm -1 appear in E||b transmittance spectra as well ||^l). We have studied the resonances at 91 and 
939 cm" 1 (E||a); 215 and 225 cm" 1 (E||b) for the samples of different thicknesses and found that while the intensity 
of 91, 939 and 225 cm -1 lines grow proportionally to the sample thickness d (that is, a = const), the intensity of the 
215 cm -1 line does not, practically, depend on the thickness (ad ~ const). Consequently, while the frequencies 91 
and 939 cm" 1 (E||a) and 225 cm" 1 (E||b) correspond to intrinsic resonances, u> = 215 cm -1 must refer to a surface 
excitation. All the observed infrared phonon frequencies together with the calculated ones are displayed in Table 1. 

NaX^Os crystals are well transparent in the frequency region between 2500 and 4500 cm -1 and below 100 cm" 1 . In 
these regions, an interference pattern was observed in (E||a) and (E||b) transmittance spectra of the samples of good 
quality. We also managed to observe the interference pattern below 100 cm" 1 in E||c transmittance of 1.3 mm thick 
sample. By measuring the distances A between the interference maxima, we have found the values of the refractive 
indexes n according to the relation 

d being the thickness of a sample. The appropriate values of e = n 2 are listed in Table 1. 

We also looked for the higher order vibrational spectra by measuring the transmittance of thick (d=0. 4-3.0 mm) 
samples in the frequency range 1000-4000 cm -1 . While no pronounced features were found in E||a and E||b polar- 
izations, sharp resonances were observed in E||c (k||a) polarization at 1930, 2858 and, possibly, 1072 and 1270 cm" 1 , 
the latter two lines being somewhat masked by the edge of a strong phonon at 955 cm -1 (see Fig. 5). 



C. Raman spectra 

Polarized room-temperature Raman spectra of NaA^Os in the spectral range 80-1000 cm" 1 are shown in Fig. 6. 
One can see immediately that the three diagonal components aa, bb, cc of the Raman scattering tensor differ markedly 
one from another which points to a considerable anisotropy of the structure. The most intense spectra were observed 
in the A g geometry a(cc)d. The intensity of the lines marked by asterisks in Bi g (i — 1, 2, 3) spectra depended strongly 
on slight variations in the sample orientation. Evidently, these lines are present due to a leakage of strong lines from 
A g geometries. We failed to assign for certain a weak feature near 100 cm" 1 in the b(ac)b spectrum overlapped by a 
strong unshifted line present in this geometry. Possibly, it comes from a leakage of a very strong line 90 cm" 1 from 
the (cc) polarization. Frequencies of the observed Raman modes together with the calculated ones are collected in 
Table 2. 

As we have already communicated pl| , besides relatively narrow lines, a broad band with a maximum near 600 cm -1 
is observed in the c(aa)c spectrum (see Fig. 6). Since this band appears under both 514.5 nm and 488 nm excitation, 
we conclude that it originates from the Raman scattering process. However, a big width of this band (213 cm" 1 ), 
practically independent of the temperature, does not permit to assign it to fundamental modes. 

We have also studied Raman spectra of Na-deficient samples Nai_a;V205 (x = 0,0.05,0.10,0.15). The most 
prominent changes occur in the A g (aa) spectrum (see Fig. 7). The 447 cm" 1 Raman line moves to higher frequencies 
with growing x. Its position shown by the vertical dashed lines in Fig. 7 is 477 cm" 1 for the sample with x = 0.15. 



4 



The maximum of the broad band moves in the opposite direction, namely, from 632 cm -1 for x = to 562 cm -1 
for x = 0.15. Such a change of the frequency difference between these two Raman bands is caused, probably, by a 
change of intermode interaction. The shape of the broad band can be approximated well by a Gaussian for all the 
values of x, its width growing from 213 cm -1 for x = to 290 cm -1 for x — 0.15. As for phonon Raman lines, 
their shape is almost Lorenzian, their width grows too. For example, the lines 177, 301 and 531 cm -1 broaden from, 
correspondingly, 11, 18 and 20 cm -1 for x = to 16, 27 and 34 cm -1 for x = 0.15. The broadening of Raman bands 
with growing x is, evidently, connected with an increase of lattice disorder. 

It is difficult to compare absolute intensities of the spectra for different x. However, certain conclusions concerning 
the relative intensities in a given spectrum can be drawn. The most prominent features are the growth of the intensity 
of the 301 cm -1 line and the appearance of a new line at 988 cm -1 for x = 0.15. All these results were obtained by 
decomposing of the observed spectrum into separate spectral contours. An example of such a decomposition is shown 
in Fig. 8. 

D. Calculations of vibrational spectra 

To obtain an information about the phonon spectrum of NaV20s throughout the Brillouin zone which is necessary 
for the analysis of the spin-phonon interaction effects, we have considered the lattice dynamics of this crystal in the 
framework of the rigid ion model. The goal of this study is to display the basic pairwise interion interactions which 
determine the main features of the measured Raman and infrared transmittance and reflectance spectra. 

From the measured large TO-LO splittings of some normal modes at the Brillouin zone center (r point), it is clear 
that long-range Coulomb forces play an essential role in the formation of the vibrational spectrum of NaA^Ov The 
potential energy of the lattice was represented by a sum of Coulomb and non-Coulomb interactions. The Coulomb 
terms in the dynamical matrix were calculated exactly by using the Ewald method. Non-Coulomb interactions 
in the form of the Born-Mayer potentials with the exponential dependence on the interionic distance r (ipij (r) = 
Aij exp(— r/pij)) were introduced between V-0 (five bonds per vanadium ion), Na-0 (eight bonds per sodium ion) 
and 0-0 neighboring ions at interionic distances shorter than 0.325 nm. Because of the nonequivalence of the oxygen 
01, 02 and 03 ions we should introduce different potentials for different types of bonding. At the initial step we 
confined ourselves by only 10 fitting parameters (instead of 34 independent force constants for the V2O5 lattice in Ref. 
p2fl ) including ion charges Z(A) (a condition of the lattice neutrality brings the relation Z(N&) + Z(01) + 2Z(V) + 
2Z(02) + 2Z(03) — about) and Aij, pij constants for V-O, Na-0 and 0-0 pairs of ions. 

The theoretical analysis of the vibrational spectra has been carried out for both proposed in the literature lattice 
structures. We did not obtain any physically well-grounded set of parameters which might provide the stable C 2v 
lattice structure. So, we discuss further in this section only vibrations of the centrosymmetric lattice with the D\\ 
space group. 

The orthogonal transformation of the atomic displacements to the symmetrized and normalized linear combinations, 
namely, 

u a (T lu ,A) = (u la (A) + u 2a (A) + u 3a (A) + u 4a (A))/2, 

u a (T lg , A) = (-u la (A) - u 2a (A) + u 3a (A) + u 4a (A))/2, 

u a (T 2u , A) = (-ui a (A) + u 2a (A) - u 3a (A) + u 4a (A))/2, 

u a (T 2gi A) = (u la (A) - u 2a (A) - u 3a {A) + u ia {A))/2, 

u a (T lu , B) = (u la (B) + u 2a {B))/V2, 

u a {T lg , B) = {-ui a {B) + u 2a {B))/V2, 
where A stands for V, 02 or 03 ions, and B — for Na or 01 ions, divides the dynamical matrix at the T-point into 
blocks corresponding to the irreducible representations of the crystal factor-group. Here, 

for a — x, Ti u — B 3u , Ti g — B 2g , T 2u = B lu , T 2g = A g , 

for a — y, Ti u — B 2u , Ti g — B 3g , T 2u = A u , T 2g = Bi g , 

for a — z, Ti u — Bi u , Ti g = A g , T 2u = B 3u , T 2g = B 2g . 

From fitting of the calculated eigenvalues of the dynamical matrix to the measured frequencies of the lattice normal 
modes, we obtained the following values of effective ionic charges: Z(V) — 2.405, Z(Na) = 0.83, -Z(Ol) = —1.22, 
Z(02) = —1.23, Z(03) = —0.98 in units of the proton charge. The apical oxygen ion 03 closest to the vanadium 
ion has the lowest charge. The average values of parameters defining the Born-Mayer potentials equal A^ = 12.2; 
0.82; 12.8 (in units of 10 3 eV), p tj =0.0181; 0.031 and 0.0212 nm for the V-O, Na-0 and 0-0 bonds, respectively. 
It should be noted that we had to correct some force constants calculated with these potentials (the relative values 
of corrections are the most significant for the distant 0-0 bonds but do not exceed 50%) to improve a comparison 
between theory and experiment. 
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The calculated frequencies of the lattice normal modes at the Brillouin zone center are presented in Tables 1 and 
2 for infrared and, correspondingly, Raman active modes. The calculated frequencies of silent A u modes equal 120, 
167 and 572 cm -1 . Acoustical properties of the lattice are defined by nine elastic constants, the predicted values of 
Cn = 17.7, C\2 = 9.7, C22 = 23.6 (in units of 10 10 N/m 2 ) are less dependent on variations of the model parameters. 

The measured components of the high frequency dielectric tensor e m remarkably differ from unity (see Table 1), 
thus a neglect of the electronic polarization is a very crude approximation in this case, in particular, when estimating 
LO-TO splittings in the T point. However, for most of the infrared active normal modes, our model gives a sat- 
isfactory description of the longitudinal macroscopic electric field induced by the vibrations of ions. A very strong 
damping of the B\ u TO mode at the frequency of 591 cm -1 may be the reason for a large difference between the 
calculated and measured LO-TO splitting (see Table 1); in the case of the B 2u TO mode at 365 cm -1 our model gives 
strongly overestimated frequency of the corresponding LO mode. Rather large discrepancies between several calcu- 
lated frequencies of the Raman active modes and the experimental data (see Table 2) clearly demonstrate that some 
significant interactions, in particular, the three body forces, greatly affecting frequencies of the bending vibrations, 
are to be included in the more thorough study of the lattice dynamics of this system. 

May be, the most interesting result of this analysis of the NaV205 lattice dynamics is the predicted soft mode 
behavior of the transverse acoustical mode at the Brillouin zone boundary (with the wave vector ko = 7r(0, 0, 1/c)) 
polarized in the ac-plane. Due to a competition of the long range Coulomb and short range non-Coulomb forces, the 
corresponding branch of the vibrational spectrum moves to the range of imaginary frequencies when approaching the 
k point thus opening a possibility to consider the NaV20s crystal as an improper virtual ferroelastic. To stabilize 
the lattice in respect to ko -excitations, we had to introduce an attractive interaction between the neighboring Vi 
and V2 (V3 and V4) ions along the a- axis with the significant bending force constant of approximately 5 N/m. The 
charge ordering in the subsystem of V ions can destroy a balance between the forces of different signs and induce a 
freezing of the soft mode atomic displacements (the unit cell doubles in the c direction, the neighboring layers shift 
in the opposite directions, in each layer the right and left legs of the vanadium ladders become nonequivalent due to 
shifts of V1-V2 and V3-V4 rungs along opposite directions in the ac-plane) as a precursor of the subsequent magnetic 
ordering with doubling of a unit cell in a-, b- and c-directions. 

IV. DISCUSSION 

A. Symmetry group of NaV^Os 

Table 3 summarizes the observed vibrational modes together with their interpretation both in centrosymmetric 
D\\ and noncentrosymmetric C\ v groups. While the former group explains naturally the experimental data provided 
one Raman and three infrared frequencies not being detected, the latter group leads to an assumption that 22 of 45 
expected Raman and 23 of 37 expected infrared modes were not detected. Moreover, only three frequencies (90 cm , 
174 and 951 cm -1 ) coincide within the experimental precision (±4 cm -1 ) in the sets of Raman and infrared modes 
corresponding to a given irreducible representation of the noncentrosymmetric group, while all the modes should be 
both Raman and infrared active in that case. We also remind once more that we failed to obtain a realistic set of 
force constants when carrying out the lattice dynamics calculations in the assumption of C 2v noncentrosymmetric 
space group. 

We consider our Raman and infrared data as well as the results of lattice dynamics calculations to support strongly 
the conclusion of the previous structural study Jl3|jl4| that the space group of NaV20s above T c = 35 K is the 
centrosymmetric rather than noncentrosymmetric CL, group. From the point of view of group it is also 
easy to explain the results of a recent 51 V-NMR study |l6| that revealed only one vanadium position at elevated 
temperatures. 

B. Atomic displacements 

As we have already mentioned in the Introduction, the structure of NaV20s looks like the structure of V2O5 
intercalated with Na. The V-0 bond lengths within the vanadium-oxygen layers are close in these two compounds 
(see Table 4). The longest bond within the layer interconnects two V2O5 units in the crystal unit cell (see Fig. lb 
where this bond is indicated as a dashed line). So, it is of a certain meaning to classify the k = crystal vibrations 
on the basis of internal vibrations of the V2O5 "molecule" (C 2v point symmetry group) split into Davydov doublets 
of the D 2 h factor group by an interaction between two "molecules" in the crystal unit cell. (A g + Bi u ), (B 2g + B 3u ), 
(Bs g + B 2u ) and (B\ g +A U ) Davydov doublets come from A\, B±, B 2 and, correspondingly, A 2 vibrations of the V2O5 
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"molecule". Splittings of these doublets can be as big as 100 cm" 1 due to Coulomb interactions and, in particular, due 
to interactions between the adjacent V ions through the common neighbors (02 ions) along the chains (see Fig. 1). 
Many of the vibrational frequencies of Na V2O5 are close to those of V2O5 p2). 

The comparison of our observed vibrational frequencies with those of V2O5 |E2] and with the results of our calcu- 
lations leads to the following assignment of the vibrational modes of NaX^Ov The V-03 stretching modes manifest 
themselves by two Davydov doublets: 951 cm -1 (B 2g ) + 939 cm" 1 (B 3u ) and 970cm" 1 (A g ) + 955 cm -1 (Bi u ). 
The Davydov splittings are relatively small in this case indicating that these vibrations associated with the strongest 
bond V-03 are really well localized. The mode frequencies are somewhat lower than the corresponding frequencies in 
V2O5 (976, 982, 994 and 975 cm -1 ), which is in accordance with longer V-03 bonds in NaV205 in comparison with 
V2O5. The following vibrations are associated with the Ol-V-03 bending modes: 177 cm -1 (A g ), 162 cm -1 (Bi u ), 
392 cm" 1 (B 2g ), 366 cm" 1 (B 3g ). 

The bridging oxygens 01 participate in the bending V-Ol-V vibrations 418 cm" 1 (B 3g ), 367 cm" 1 (B 2u ), 
447 cm -1 (A g ), 468 cm -1 (B\ u ). The stretching V-Ol-V vibration is placed at 420 cm" 1 (A g ) and involves mainly 
the motion of the V atoms along the a-axis. 

The modes at 683 cm -1 (B% g and B\ g ) and 582 cm" 1 (B 2u ) correspond to the stretching V-02 vibrations along 
the 6-axis while those at 550 cm -1 (B 2g ), 533 cm" 1 (A g ) and 526 cm -1 (B 3u ) — to the bending ones. 

Most of the remaining modes can be described in terms of external modes of the V2O5 units. Thus, the modes at 
186 cm -1 (B2 g ), 169 cm -1 (i?3 S ) and 90 cm -1 (A g ), correspond to the relative translations of the two V2O5 units 
within the crystal unit cell along the a-, b- and, correspondingly, c-axes. As the V2O5 units are bound together along 
the 6-axis, these modes can be considered as relative translations of the neighboring (V05)„ chains. The B ig mode 
at 174 cm" 1 (03 ions move along the &-axis) and the A g mode at 301 cm" 1 (02 ions move along the c-axis) correlate 
with in-plane and, correspondingly, out-of-plane chain bending vibrations. The B 2g mode at 141 cm -1 and B 3u mode 
at 91 cm -1 are connected with the rotation of the chains around the &-axis. 

The following modes involve mainly the displacements of Na atoms: 225 cm" 1 (B 2u ), 251 cm" 1 (B 3u ), 
181 cm" 1 (Biu). 

C. Spectra of electron excitations 

With obtained values of the effective charges, we estimated the crystal field energies of the 3d electron localized at the 
V 4+ ion site. The crystal field parameters B% =1360+2090G, B\ =2020-1590G, B\ =820+640G, B\ =610+1430G, 
B\ =-1810-3690G, B\ =33+144G, B\ =3650+8180G, B\ =4070+7420G cm" 1 for the Vi, V 3 sites were calculated in 
the framework of the exchange charge model j27J (for the V2, V4 sites B 1 and B^ parameters change signs, the first 
terms correspond to point charge contributions, the Stevens normalization is used). The scaling factor G determines 
the strength of the exchange charge field. We estimated this phenomenological parameter of the model (G = 4) by 
fitting the total crystal field splitting to the width of the V-3d bands presented in Ref. Jig] . In this case, the effective 
crystal field provides the following energy level pattern of the V 4+ ion : (A"), 1.10 (A'), 1.18 (A"), 3.39 (A 1 ) and 
4.78 (A 1 ) eV (irreducible representations of the G s point group, corresponding to the space symmetry of the electron 
wave function, are given in brackets, additional shifts of the crystal field levels due to the spin-orbit interaction and 
the electrostatic field of a hole at the neighboring vanadium site are less than 0.025 eV). The ground state wave 
function is the d xy orbital with small admixture of the d yz orbital as it has been already pointed earlier fl3|| , and the 
sequence of the excited states is in agreement with results of the band structure calculations as well [|l3| . 

The strong absorption of light (k||c) with the E||a as well as with E||b polarizations was observed in the region of 
1.0-1.2 eV plj . Both magnetic dipole and induced electric dipole d — d transitions in the odd crystal field are allowed 
between the A" states for E||a, and between the A" and A 1 states for E||b. Thus, in accordance with the results on 
the crystal field energies given above, the observed broad optical bands in the region of 1 eV can be interpreted as 
phonon assisted d — d transitions without any additional suppositions about the broken symmetry between legs of 
vanadium ladders p8| . 

The next step towards the detailed description of the spectra of electron excitations is to construct molecular 
orbitals for the [V20] 7+ "molecule" having G2^ point symmetry, from the vanadium d-orbitals and oxygen p-orbitals. 
The vanadium ground state wave function d xy yields the non-bonding a 2 orbital as well as bonding and antibonding 
62 molecular orbitals, namely, 
a 2 : [d X y(L) + d xy (R)} 
b 2 : [d xy (L) - d xy (R)} + p y 
b* 2 : \d xy (L) - d xy (R)]-p y . 

Here, a 2 , b 2 denote irreducible representations of the C 2v point group, the asterisk refers to an antibonding orbital, 
L and R stand for the vanadium sites at the left and, correspondingly, right sides of a ladder rung. The highest 
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filled orbital being a^, the 0,2 — > 62 electronic transition allowed in E||a polarization could account for the low- 
frequency absorption band observed just in this polarization. Quantum-chemical calculations are necessary to verify 
this qualitative interpretation. 

D. Fano resonances with a continuum 

The asymmetric lineshapes of the infrared active modes at 91, 150, 939 cm -1 in E| |a polarization bring into evidence 
a strong interference between these modes and a continuum observed just in this polarization. This interpretation is 
supported by the fact that the spectral line near 91 cm -1 becomes perfectly symmetric when the continuum absorption 
vanishes in this spectral range below the phase transition temperature T c = 35 K. In our earlier work || we argued 
that these changes are connected with the opening of a gap in the magnetic excitation spectrum at T c , the observed 
continuum being due to two-magnon absorption. 

However, such a straightforward interpretation is no more valid in the case of the space group D^l- It should be 
revised, taking into account possible electronic excitations in this frequency range, as discussed in the previous section, 
and a charge ordering at the transition temperature. 

E. Higher order infrared vibrational spectra 

Two- and three-phonon absorption occurs due to anharmonicity of crystal vibrations. It is continuous, displaying 
singularities corresponding to critical points of the Brillouin zone. Leaving the detailed analysis of multiphonon bands 
to a future publication, we discuss here only sharp lines observed in E||c absorbance spectrum (Fig. 5). They are 
listed in Table 5 together with their tentative assignment, using symmetry allowed combinations of T-point phonons 
observed in our first-order spectra. The coincidence of the observed and combinational frequencies is within the 
precision of our measurements. 

The strongest narrow peak at 1930 cm -1 corresponds, according to this assignment, to a sum of the components 
of the Davydov doublet originating from the V-03 stretching vibration. This stretching mode is well localized which 
results in its small dispersion over the Brillouin zone, thus delivering a narrow two-phonon band, in accordance with 
the experimental observation. 

V. SUMMARY 

We have performed a thorough spectroscopic study of far infrared reflection and transmission as well as Raman 
scattering of a'-NaX^Os single crystals in the high temperature phase (above T c = 35 K). Far infrared spectra were 
obtained for E||a, E||b and E||c polarizations of incident light. Diagonal (aa), (bb), (cc) as well as nondiagonal (ab), 
{be), (ac) components of the Raman scattering tensor were investigated. We report five infrared active modes in 
E||a polarization, four — in E||b polarization and six — in E||c polarization. Eight Raman active modes have been 
detected for parallel polarizations of incident and scattered light (aa). (bb), (cc). The (ab), (ac) and (be) Raman 
geometries delivered three, seven and, correspondingly, five modes. These results are in a much better agreement with 
the recently proposed centrosymmetric space group DtA (Pmmn) for the high-temperature phase of NaV20s than 
with earlier accepted noncentrosymmetric space group C\ v (Pmn2i). We have also performed the lattice dynamics 
calculations, in the framework of the rigid ion model, for both proposed in the literature structures of NaV2 05. 
We failed to obtain any physically well-grounded set of parameters providing a stable C\ v lattice structure. Thus, 
our infrared and Raman experimental data, as well as the results of lattice dynamics calculations support strongly 
the conclusion of the previous structural study Jl3|,|l4| that the space group of NaV20s above T c — 35 K is the 
centrosymmetric D^f. rather than noncentrosymmetric C\ v group. 

This conclusion leads to important physical consequences. In particular, it requires to revise the interpretation 
of one-dimensional magnetic properties of NaV20s and of the phase transition at 35 K considered earlier as an 
ordinary spin-Peierls transition. An interpretation of the earlier observed broad bands in the near and far infrared 
absorption J2l],|| needs to be reconsidered too. 

Using the values of effective charges obtained in the process of lattice dynamics calculations and fitting the total 
crystal field splitting to the width of the V-3<i bands |l|] , we estimated the crystal field energies of the 3d electron 
localized at the vanadium site. It follows from this estimate that the observed near infrared broad band absorption 
of NaV2C>5 |2lJ can be interpreted as phonon assisted d — d transitions. We speculate that the far infrared E||a 
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polarized absorption continuum might be associated with electron excitations of [V2 0] 7+ rungs in the crystal field of 
Ci v symmetry. 

Strongly asymmetric spectral lines observed in E||a absorbance spectra of NaV20s bring into evidence a strong 
interference between relatively narrow phonon lines and the underlying continuum. This points to an interaction 
between crystal vibrations and magnetic or electronic excitations. The detailed physical interpretation of the observed 
phenomenon depends on the nature of the far infrared E||a polarized continuum which requires a special investigation. 

In conclusion, we reported also some preliminary results on higher order vibrational spectra of NaV2 0s appearing 
due to anharmonicity of lattice vibrations. 
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FIGURE CAPTIONS 



FIG. 1. The structure of NaV^Og. (a) A stereometric projection. Oxygen and vanadium atoms are at the corners of 
and, correspondingly, inside the pyramids. Sodium atoms are presented as balls, (b) ab projection. Apical oxygen 03 atoms 
(situated above or below the corresponding V atoms) are not shown. A dashed line indicates the longest V-02 bond (0.199 nm). 
a&-projection of the crystal unit cell is shown by a thin solid line. 

FIG. 2. Room temperature far-infrared reflectivity spectra of NaV20s. Open circles represent experimental data. Solid lines 
are the result of the fitting (see the text). 

FIG. 3. Absorption coefficient a in the region of low frequency absorption bands at room temperature. An arrow indicates 
a Fano-type resonance. It is shown separately in the inset. Open circles represent experimental data. Solid line in the main 
figure was calculated using the parameters obtained from the fitting of the reflectance spectrum. Solid line in the inset is a 
result of the fitting by the expression (5) with shown as a dashed line. 

FIG. 4. Fano resonance near 91 cm -1 at the temperature of 40 K (open circles) and its fitting by the equation (5) with 
parameters uj r = 90.7 cm -1 , 7 = 0.2 cm -1 , q = —1.0, 0:3(^2) = 270 cm -1 , olq/olb = 0.3 (solid line). The temperature 
dependence of the Fano parameter q is given in the left inset. The right inset presents the absorption spectrum in the 
vicinity of 91 cm -1 line at 6 K taken with the resolution of 0.05 cm -1 (open circles) and its fitting by the Lorenzian with a 
FWHM=0.10 cm- 1 . 



FIG. 5. Absorbance spectrum of NaV20s in the region of multiphonon bands at room temperature. 



FIG. 6. Room temperature Raman spectra of NaV20s. Asterisks mark A g lines seen in Bi g (i = 1, 2, 3) spectra. 



FIG. 7. Room temperature Raman spectra of Nai-^V^Os for different x. 

FIG. 8. Decomposition of the room temperature c(aa)c Raman spectrum of NaV20s into separate contours (dashed lines). 
The sum of these contours shown by a solid line approximates the experimental spectrum (circles) rather well. 
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TABLE I. Infrared active vibrational modes and dielectric constants of NaV20s (all frequencies are in cm x ) 



Polarization, 



Observed 



mod© 


Transmission 


Reflection 






symmetry 


T = 


40 K 




T = 


300 K 








T 


= 300 K 














LOTO 


fTO 


LOTO 


7to 




a 


LOTO 


7to 


LOhO 


JLO J 


• 10 3 






LOTO 


IjJlO 


E||c 












7.5 


162 


5.2 


165 


5.7 


45 


3.9 


7.7 


216 


219 


Bi u (Pmmn) 












±0.2 


179 


8.4 


212 


8.3 


130 






232 


256 


or 




























298 


298 


B 2 (Pmn2 1 ) 
















oo.U 


/ICQ 
4oo 


oo.U 


ZO 






430 


430 
















591 


119.9 


594 


119.0 


6.8 






ooy 


con 
















760 


59.4 


762 


59.3 


1.6 






691 


716 
















955 


2.5 


1017 


3.0 


39 






961 


1036 


E||b 


178 


4 


175 


12 


5.2 


10.2 


175 


8.3 


180 


8.4 


39 


4.9 


9.5 


141 


173 


B2u(Pmmn) 


225 


1 







±0.2 


±0.2 
















240 


266 


or 


371 




367 


16 






365 


12.8 


378 


13.3 


60 






388 


483 


B 1 (Pmn2 1 ) 


594 


13 


582 








584 


29.5 


769 


29.0 


271 






578 


747 


E||a 


91 c 








9.6 


15.0 












7.7 


15.8 


111 


126 


B3 U (Pmmn) 


140 b 




145 b 




±0.3 


±0.6 


153 


33.0 


154 


34.1 


46 






130 


177 


or 


254 




251 








251 


7.3 


252 


8.2 


9 






227 


276 


Ai(Pmn2i) 




























493 


534 




531 


18 


526 


53 






525 


39.5 


621 


52.4 


208 






538 


653 






























742 


808 








939 d 






















955 


957 



Calculated 
(Pmmn) 



a o;i = 3200 cm" 1 , LU2 = 40 cm" 1 . 
b Asymmetric line. 

90.7 cm" 
= 939 cm" 



c Fano-typc resonance: lo t 
d Fano-type resonance: uj r 



, 7 = 0.2 cm 1 , q = —1.0, ao/as = 0.3. 
7 = 1.0 cm" 1 , q — 1.1, Qo/qs = 0.2. 



TABLE II. Room temperature Raman frequencies (cm x ) for NaV20s 



aa, bb, cc 


ab 




ac 




6c 


A\ g (Pmmn) 


Big(Pmmn) 




B2g(Pmmn) 


B 3 g 


(Pmmn) 


or A 1 (Pmn2 1 ) 


or Bi(Pmn2i) 


or B 2 (Pmn2i) 


or A 2 (Pmn2i) 


Observed Calc. (Pmmn) 


Observed Calc. (Pmmn) 


Observed Calc. (Pmmn) 


Observed 


Calc. (Pmmn) 


90 (cc, aa) 91 


174 191 


141 


129 


169 


149 


177 (aa, bb, cc) 226 


295 288 


186 


193 


257 


239 


233 (bb) 319 


683 679 


225 


296 


366 


262 


301 (aa,bb) 362 




392 


332 


418 


396 


420 (66, cc) 439 




429 


410 


683 


685 


447 (aa) 518 




550 


619 






533 (aa, bb) 626 






798 






970 (cc, aa) 964 




951 


961 
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TABLE III. Comparison of experimentally observed Raman and infrared modes with the expected ones within centrosym- 
metric D 2 jJ and noncentrosymmetric C 2 „ space groups (mode frequencies are in cm -1 ) 

Pmmn(DW) Observed modes Pmn2\(C\ v ) 



8A g (aa y bb,cc) 90 177 233 301 420 447 533 970 \ 

15Ai (aa, 66, cc;E||a) 

7B 3t ,(E||a) 91 a 145 251 526 939 / 



W lg (ab) 
4B 2tl (E||b) 


174 
175 


295 
225 a 


683 
367 


582 








\ 
/ 


7Bi (afe;Ej|b) 


8B 2g (ac) 
7B Ul (E||c) 


141 
162 


186 
179 


225 
468 


392 
591 


429 
760 


550 
955 


951 


\ 
/ 


15B 2 (ac; E]|c) 


5B 3g (6c) 
3A U 


169 


257 


366 


418 


683 






\ 
/ 


8A 2 {be) 



a Observed below 200 K. 



TABLE IV. The bond lengths (nm) in NaV 2 5 and V2O5 



Bond 



NaV 2 5 [IS] 



V 2 5 [22] 



V-03 

V-Ol 

(V-02)x2 

V-02' 

V-03' 



0.161 
0.183 
0.192 
0.199 
0.318 



0.158 
0.177 
0.188 
0.202 
0.278 



TABLE V. Multiphonon bands in NaV 2 Os observed in E||c polarization 



Observed bands, cm 1 


Combination of phonons, cm 1 


1072 (B la ) 


550 (B 2g ) 


+ 526 (B 3u ) 


= 1076 (Biu) 


1270 (Sit,) 


683 (B 3g ) 


+ 582 (B 2 „) 


= 1265 (BiJ) 


1930 (B lu ) 


970 (A g ) 


+ 955 (B lu ) 


= 1925 (Bi«) 


2858 (BiJ) 


3 x 955 (Biu) 




= 2865 (B lu ) 


2901 (Bi u ) 


2 x 970 (A g ) 


+ 955 (Biu) 


= 2895 (BiJ) 
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